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Abstract 

The equivalence of thermodynamic results in the canonical and the 
microcanonical ensembles has been questioned in some calculations for 
spin models with long-range interactions. We show that these claims 
of inequivalence are related to an inadequate definition of the indepen¬ 
dent (density) variables in the microcanonical ensemble. We illustrate 
this point with the example of a simple spin-1 ideal paramagnet, and 
then revisit the original calculations of Barre, Mukamel, and Ruffo, 
for a mean-field spin-1 Blume-Capel model. If the microcanonical en¬ 
semble is defined in terms of adequate density variables, we show that 
there is no disagreement with the calculations in the usual canonical 
ensemble (with fixed thermodynamic field variables). 

The equivalence of ensembles, which is one of the hallmarks of modern 
statistical physics urn has been challenged by Barre, Mukamel and Ruffo 
[3], in the calculations for a simple mean-field, fully connected, spin-1 Hamil¬ 
tonian. This work motivated a number of remarks on similar thermodynamic 
discrepancies between results in the usual canonical and microcanonical en¬ 
sembles associated with spin model systems with long-range interactions [1]. 

We show that this claimed inequivalence of ensembles is due to an inade¬ 
quate definition of the independent variables in the microcanonical ensemble. 
A proper definition of the microcanonical ensemble is illustrated by analytic 
calculations for a spin-1 ideal paramagnet. We then revisit the original work 
of Barre and collaborators for a mean-field version of the spin-1 Blume-Capel 
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model, which is known to display a temperature-anisotropy phase diagram 
with second and first-order transition lines and a tricritical point. If the mi- 
crocanonical ensemble is defined in terms of the adequate density variables, 
we show that there is no disagreement with the well-known solutions in the 
usual canonical ensemble (with fixed temperature and thermodynamic held 
variables). According to Griffiths and Wheeler [5], we distinguish between 
“density” variables, given by the ratio of thermodynamic extensive quanti¬ 
ties, and “held” variables, which “take on identical values in two phases that 
are in thermodynamic equilibrium with each other”. 


1 Ideal paramagnet of spin-1 

In this Section, we use a spin-1 ideal paramagnet to give a very simple exam¬ 
ple of the choice of variables in the usual canonical and the microcanonical 
formalisms. Consider a spin Hamiltonian given by 

N N 

n = -Hj2S) + D Y, s h (i) 

3 = 1 3=1 

where H and D are “held variables”, and Sj = —1, 0, +1, for j = 1, 2,..., N. 


1.1 Calculations in the canonical ensemble 


In the canonical ensemble (with hxed values of the held variables, H , D , and 
temperature T), we have the partition function 

r n n i 


Z = Z(T, H, D ) 


exp pHj2Sj-P D J2 s i 

{Sj} l 3=1 3=1 - 


= [1 + 2 exp (—/ 3D ) cosh (f3H)] N , (2) 

where /3 = 1 jk^T. The associated thermodynamic potential per site is given 
by 

g = g(T , H, D) = ~ In [1 + 2 exp (-0D) cosh (/ 3H)\ . (3) 


This is a fundamental equation E m the 'T is j. ^oviiucuji'jir. 

this thermodynamic potential, we obtain the equations of state, 

(rp tt t -)\ (dg \ sinli (0H) 

m = m(T,H,D) = - I — ' 


T,D 


l exp 


(0D) + cosh (0H) ’ 


(4) 
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where the thermodynamic density m is the dimensionless magnetization per 
site, m = (y2f =1 Sj \ /N, and 


q = q(T,H,D) 


dg \ cosh (/3H) 

dD ) T H | exp (/3D) + cosh (j3H) ’ 


(5) 


where the density q is a dimensionless “quadrupole”, q = /N- 

The entropy per site is given by the remaining equation of state, 


s = s(T,H,D) = -(^) = k B In [1 + 2 exp (-/3-D) cosh {/3H)\ + 


H,D 


k 




/3D cosh (/ BH) 


pH sinh (pH) 


2 exp (PD) + cosh ( pH ) \ exp (PD) + cosh (PH) 


( 6 ) 


In this usual canonical representation, the thermodynamic densities m, q, 
and s, are obtained as functions of the fields T, H, and D. The meaning 
of D as a thermodynamic field can be fully appreciated if we use the spin-1 
model to mimic a binary lattice gas of particles, in which case D is related 
to a chemical potential. 

It is now interesting to analyze some features of this ideal spin-1 param- 
agnet. Using equations (j4j) and (J5]) , the entropy can be written as 


1 D H 
s = ~j9 T —q — m 


(7) 


which shows that, in this schematic model, the held terms in the effective 
Hamiltonian do not contribute to the “magnetic internal energy”. In other 
words, the average value of the spin Hamiltonian is a kind of “magnetic 
enthalpy”, which cannot be taken as the internal energy |6]. Equation (|7D 
is just the analogous of the Euler relation for a simple fluid system. It is 
then straightforward to use this Euler relation, and the expression for the 
free energy, g = g (' T , H , D), to write the differential form 


ds 


Y d q - -f dm ’ 


( 8 ) 


which shows that the entropy per particle s can also be written as a function 
of the densities m and q. In this special case, we can perform the appropriate 
calculations to write the analytic expression 

s = s (q,m) = —ks (1 — q) In (1 — q) — ks~ (q — m) In (q — m) — 
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- k B ^{q +m) In (q + m) + k B qln2, (9) 

which is equivalent to a fundamental equation in the microcanonical ensemble 
(in terms of the densities m and q, and the number of sites N ; note the fixed 
value, u = 0, of the internal magnetic energy in this simple paramagnetic 
model). Incidentally, this expression has the same form as the “entropy to be 
maximized”, given by equation (5) in the article of Barre and collaborators 

0 . 


1.2 Calculations in the microcanonical ensemble 


We now show that the same expression for the entropy per particle, s = 
s(q,m ), given by equation (jUJ), can also be obtained in the context of a 
properly defined microcanonical ensemble, in terms of the densities m and q. 

In this simple problem, the number of microscopic configurations is given 
by 

AT! 

n = n (n + , n 0 , n_) = —-—-—(io) 

V + ’ ; N + \N 0 \NJ. V ; 

where N is the total number of spins, and N +1 Nq, and 1V_, are the numbers 
of spins up, zero, and down, respectively. Taking into account that 

N = N + + N 0 + 7V_, (11) 


and that, 

m = 1 (AT + - JV_), q= l(N + + N_), (12) 

we can write as a function of N, m, and q. In the thermodynamic limit, 
we show that 

s = s(q,m)= lim —ks In 12, (13) 

N —>oo N 

where s = s (q, m) is a fundamental equation in the entropy representation, 
given by the same expression, equation (I9l) , which we have already obtained 
in the context of the canonical ensemble. From this form of the entropy, we 
obtain the equations of state in the entropy representation, 


D / ds\ H (ds \ 

T \ d( l) m' T 


(14) 


from which we regain the usual expressions of m and q in the canonical 
ensemble, given by equations (J4| and (J5]). 
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At this point we have shown the full equivalence between the canonical 
and the microcanonical formulations of this simple problem. We can also re¬ 
sort to an alternative (and more general) formulation in the microcanonical 
ensemble, which turns out to be useful for dealing with short-range interact¬ 
ing systems as well. 

In the alternative formulation, we rewrite the canonical partition function 
in the form 


Z = Z(T,H,D) = J2 f dm [d q s( m - i £ S, ) s( q - i x 

{s j y J J V j=i ) \ j= i / 

x exp (/ QHNm — (5DNq ) , (15) 

which can also be written as 

Z = J dm J dqQ(q,m,N) exp(/3HNm — f3DNq) , (16) 


where hi (q, m, N ) is the number of microstates with fixed values of q, m, 
and N. We now introduce convenient integral representations of the delta 
functions. In the thermodynamic limit, we have 

~ exp [Nf (q, m- k x , k 2 )\, (17) 


with 


/ = kyn + k 2 q + In [1 + 2 exp (— k 2 ) cosh (£q)] , 


(18) 


where the parameters k\ and k 2 come from the saddle-point equations, 


df_ 

dk\ 



(19) 


The microcanonical entropy is given by s = s(q,m) = ksf, with k\ and 
k 2 obtained from the extremization of the function / ( q,m ; ki, k 2 ). It is not 
difficult to show that k\ = —/ 3H , and k 2 = /3D } so that equation (IT8l) leads 
to the same expression of equation ([6]). Again, we recover the known results 
in the canonical ensemble. 
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2 The spin-1 Blume-Capel model 


We now consider the special spin system analyzed by Barre and collaborators 
[3], given by the fully connected spin-1 Hamiltonian 


n 


_j_ 

2 N 



N 

+ DY.SI 


3 = 1 


( 20 ) 


where J and D are positive parameters, and Sj = +1, 0, —1 for all sites. In 
the canonical ensemble (fixed values of the held variables T and D), we have 


Z = Z(T, D, N ) 


E-p 

{Sj} 


PJ 
2 N 



N 


3 =1 


( 21 ) 


Using a Gaussian identity, we write 


Z = 



'‘Too 


dy exp [—(91V g (T, D; 


( 22 ) 


where 

9 (T, D ; y) = ^y 2 - ^ In [1 + 2 exp (-pD) cosh (pjy)\. (23) 

In the thermodynamic limit, the free energy g = g (T, D ) comes from the 
minimization of (T, y) with respect to y. We then have 


g(T,D) = g(T,D-y), 


(24) 


with 


y = 


sinh (pjy) 


- exp (PD) + cosh ( pjy) ’ 


(25) 


and the proviso about the existence of multiple solutions for y as a function 
of the thermodynamic fields T and D. 

We now expand g (T,D\y), given by equation f[23|) . as a power series of 
y, with field-dependent coefficients. From this Landau-like expansion, we 
obtain the well-known results for the critical line in the D — T space, 


D = ksT In 


2 J 

k^T 


- 2 


(26) 
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with 1/3 < ksT/ J < 1, and a tricritical point at k^T / J = 1/3 (and D/ J = 
(In4) /3). For ksT / J < 1/3, there is a line of first-order transitions, with the 
coexistence of a paramagnetic and an ordered phase at well-defined values of 
the thermodynamic fields T and D. As it is usual in the mean-held models, 
the location of this first-order border conies from a proper application of a 
Maxwell construction. 

Some simple manipulations lead to the usual equations of state, for the 
entropy and the quadrupole moment, in the T — D representation, 

f dg\ cosh (/ BJy) 

q \dD J T i exp (/3D) + cosh {fBJy) ’ 

and 

S = -(f) = k B In [1 + 2 exp (—(3D) cosh (/3Jy)\ + — -^y 2 , (28) 

where y comes form the solutions of equation (j25j) . It easy to see that y 
corresponds to the magnetization m per site, Sj) /N, but we have to be 
careful in the region of multiple solutions (and first-order transitions). 


2.1 Calculations for the Blume-Capel model in the mi- 
crocanonical ensemble 

As in the work of Barre, Mukamel, and Ruffo [3], we initially write the 
number of microstates, 

n = n(N + ,N 0 ,N_) = (29) 


where N = N + + A^o + A r _ is the fixed number of spins, 
quadrupole density, 

N + + N 


We now £x the 
(30) 


and the internal magnetic energy per site, 


u = 



N_) 2 


J f N + — N_ \ 2 

2 V ) 


(31) 
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and write f2 as a function of N, u, and q. In the thermodynamic limit, the 
entropy is given by 

s — s(u,q)— lim — fc^lnfh (32) 

IV—>oo N 


In this particular mean-held model, the internal energy u is associated 
with the magnetization m, 


( N + -N \ 2 

2 V N ) 



(33) 


so that fixed values of u also correspond to fixed values of m 2 . Therefore, 
the entropy s = s (u, q) is given by 

s — s (u, q) — —k B (1 — q) In (1 — q) — k B - [q — m] In [q — m] — 

- k B ^ [q + m\ In [q + m] + k B q In 2, (34) 

which is identical to equation ()9|) , for the ideal spin-1 paramagnet, but we now 
have to take into account that m is a function of u, m = m(u), according 
to the dehnition of the internal energy of this mean-held model, given by 
equation (j33lh Note that u < 0. Also, note the symmetry with respect to 
±m. 

Given the fundamental equation for this system, s = s(u,q), we write 
the equations of state in the entropy representation, 



It is straightforward to show that these equations of state lead to the same 
expressions for the internal energy u (related to the magnetization m) and 
the quadrupole moment q, as a function of T and D, which we have already 
obtained in the context of the usual canonical ensemble. Numerical calcu¬ 
lations for the hrst-order transition border just confirm these findings. In 
the region of multiple solutions of the mean-held equations, Maxwell's con¬ 
struction should be properly used to choose the minima of the canonical free 
energy in terms of D and T (which is equivalent to recovering the convexity 
of the entropy in the microcanonical ensemble). In contrast to the work of 
Barre, Mukamcl, and Ruffo [3], these calculations indicate that there is full 
equivalence of ensembles. 



2.2 Alternative formulation of the microcanonical en¬ 
semble 

The equivalence of ensembles can also be shown if we resort to an alterna¬ 
tive, and more general formulation of the microcanonical ensemble. In this 
formulation, we introduce two delta functions to fix the energy u and the 
quadrupole density q, and rewrite the canonical partition function of the 
Blume-Capel model, 



x exp \/3Nu — /3DNq] = J du j dqQ. (■ u , q, TV) exp [— j3Nu — f3DNq ], (36) 

so that Q (u, q, N ) is the number of microstates with fixed values of u , q, and 
N. 

Using an integral representation for the delta functions, it easy to write 


+200 +200 


(u, q, N ) = 


dki f dk‘ 


J 2ni J 2ni 
—200 —200 


: exp [Nuki + Nqk 2 ] x 


x 


k\JN\ 1/2 
2t r J 


+oo 


dy exp { — -kiJNy 2 + N In [1 + 2 exp (—k 2 ) cosh (k\Jy)] }. 


In the thermodynamic limit, we have 

fl (u, q, N ) ~ exp [Nf (u, q; y, h, k 2 )] , 

with 


(37) 

(38) 


/ (u, q; y, ki, k 2 ) = kiU+k 2 q--kiJy 2 +\n [1 + 2 exp (~k 2 ) cosh (faJy )\, (39) 


where the parameters k \, k 2 and y come from the saddle-point equations, 

(40) 


df df = df = § 


dk\ dk 2 dy 
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The microcanonical entropy is given by s = s (u, q) = ksf, from which we 
write the equations of state in the entropy representation, 1/T = ( ds/du) q 
and D/T = (ds/dq) u . Using these equations, we have 


J _^2 ~ = sinh (pjy) 

2 ^ | exp (/ 3D ) + cosh (/ 8Jy) ’ 


and 


cosh (/3Jy) 

^ | exp (/3D) + cosh (f3Jy) ’ 

which are the known expressions in the canonical ensemble. 


(41) 


(42) 


3 Conclusions 

In conclusion, we have shown that claims of inequivalence of ensembles in a 
class of mean-held spin models can be attributed to an inadequate definition 
of the independent (density) variables in the microcanonical formalism. As 
an illustration of an adequate formulation, we performed some analytic cal¬ 
culations for a spin-1 ideal paramagnet. We then revisited the calculations 
of Barre, Mukamel, and Ruffo [3]. If the microcanonical ensemble is defined 
in terms of the adequate density variables, we have shown that there is no 
disagreement with the well-known solutions in the usual canonical ensemble 
(with fixed thermodynamic held variables). 
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